Dzyaloshinsky-Moriya Anisotropy in the Spin-1/2 Kagome Compound ZnCu3(OH)eCl2 
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We report the determination of the Dzyaloshinsky-Moriya interaction, the dominant magnetic 
anisotropy term in the kagome spin-1/2 compound ZnCu3(OH)eCl2. Based on the analysis of the 
high-temperature electron spin resonance (ESR) spectra, we find its main component \D Z \ — 15(1) K 
to be perpendicular to the kagome planes. Through the temperature dependent ESR line-width we 
observe a building up of nearest-neighbor spin-spin correlations below ~150 K. 

PACS numbers: 71.27.+a, 75.30. Gw, 76.30.-v 



Among all geometrically frustrated networks, the spin- 
1/2 corner-sharing kagome antiferromagnetic lattice has 
been at the forefront of the quest for novel quantum phe- 
nomena for the past two decades [lj]. Various competing 
ground states with minute energy difference have been 
proposed 0, 0] and are still under active consideration. 
The understanding of the ground state and of the low- 
lying excitations therefore still appears as a pending the- 
oretical issue, even for the simplest nearest-neighbor (nn) 
isotropic Heisenberg case. In addition, minor deviations 
from this model, often met in experimental realizations, 
such as magnetic anisotropy 0, 13] and spinless defects 
[H], may crucially affect the fundamental properties of 
the kagome antiferromagnets. 

In the pursuit of finding a suitable compound which 
would reflect intrinsic kagome properties, the min- 
eral herbertsmithite, ZnCu3(OH)gCl2, has recently been 
highlighted 0| as the first "structurally perfect" realiza- 
tion of the spin-1/2 kagome lattice [inset (b) to Fig. Q]. 
It features Cu 2+ -based kagome planes, separated by non- 
magnetic Zn 2+ . It lacks any long-range magnetic order 
(LRO) down to at least 50 mK [7| despite sizable nn ex- 
change J w 190 K Its bulk magnetic susceptibility 
Xb monotonically increases with decreasing T H fioj] . at 
variance with numerical calculations 0, [U, [lj It was 
argued that this increase is due to magnetic anisotropy of 
the Dzyaloshinsky-Moriya (DM) type fli |. Dy ■ S { x Sj, 
and/or due to impurities [9j, llOj. Various experiments 
indeed evidenced 4 — 7% of Cu 2+ /Zn 2+ inter-site disor- 
der in the samples available at present 0, Ql], Qjl, 0]. 
This explains why Xb differs significantly from the local 
susceptibility of the kagome planes Xk measured by NMR 
ljl . lq . The latter exhibits a pronounced decrease below 
~ J/2 and shows a finite plateau at low T . Together 
with ungapped excitations 3, IE 3] this points to the 
absence of the singlet-triplet gap in this compound. 

The absence of the spin gap leads to the important 
question whether this is an intrinsic feature of the kagome 



physics or is it related to additional terms in the Hamil- 
tonian of ZnCu3(OH) 6 Cl2. The presence and role of 
the DM interaction, highlighted by Rigol and Singh 
have been questioned on several occasions 0, Ql| 
Such perturbation to the isotropic Hamiltonian can in- 
deed drastically affect the low-T properties by mixing dif- 
ferent states and can even stabilize LRO, as in jarosites 
0, EH- In order to be able to inspect this crucial issue it 
is of utmost importance to experimentally determine its 
magnitude in ZnCu3(OH)gCl2. 

Using electron spin resonance (ESR), in this Letter we 
provide evidence of a sizable magnetic anisotropy, which 
we argue to be the DM interaction. In addition, through 
the temperature dependence of the line-width we show 
that spin correlations build up below ~150 K and relate 
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FIG. 1: (color online). ESR spectra of Zn^C^-^OH^Cb, 
measured at 326.4 GHz, and their fits to the powder-averaged 
Lorentzian (solid lines). The dashed line gives an axial powder 
spectrum for the same (/-factor anisotropy (g\\ = 2.25, g± = 
2.14) in the case of small magnetic-anisotropy broadening. 
Inset (b): DM vectors connecting the nn Cu 2+ cites (orange) 
through OH~ groups (blue, gray), with D z and D p as the 
out-of-plane and the in-plane components. 
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the low-T behavior to inter-site mixing defects. 

Determined by the imaginary part of the dynami- 
cal susceptibility x"( w )j ESR spectra, I(ui) oc %"(q — > 
Q,u>) oc X^di (S' + (i)5-(0))exp Mt / r , reflect the T- 
evolution of the spin correlation function (5 + (i)S _ (0)) 
in the direction perpendicular to the applied magnetic 
field H. Kramers-Kronig relations demonstrate that the 
integrated ESR intensity xesr gives a locally measured 
static susceptibility. Further, a finite ESR line-width pro- 
vides a direct measure of a finite magnetic anisotropy. 

In Fig. Q] we show typical derivative ESR spectra of 
polycrystalline ZnxCu4- x (OH)6Cl2 (x — 0.5 — 1) synthe- 
sized as in Ref. |6J], which were recorded at 326.4 GHz. 
Using a reference sample, in ZnCu3(OH)gCl2 we find 
Xesr = 1-0(2) emu/molCu at room temperature (RT), 
which is in good agreement with the bulk susceptibil- 
ity Xb = 1-1 emu/molCu. This proves that the Cu 2+ 
spins on the kagome lattice are detected by ESR 
The T-dependence of xesr is presented in Fig. [2ja). Its 
monotonic increase with decreasing T closely resembles 
Xb down to 5 K and differs from the peaked kagome 
planes susceptibility Xfe This demonstrates that 

ESR detects both the Cu 2+ spins on the kagome lattice 
and those on the inter-layer Zn sites, the latter result- 
ing from the Zn 2+ /Cu 2+ inter-site disorder. Despite the 
two detected Cu 2+ sites with likely different g-factors 
($9 "$ 0.2), only a single ESR line is detected. This is 
expected [19( because the exchange coupling J' ~ 10 K 
between the two Cu 2+ species [2Jj] is larger than their 
difference in Zeeman energy, Sg^sH < 1 K. 

In order to fit the spectra appropriately, one has to 
consider both the g-factor anisotropy and the line-width 
broadening by the magnetic anisotropy. The former is 
given by g(6) = ( 3 2 cos 2 (0 - B g ) + gl^{9 - g )f/\ 
where 9 is the angle between H and the normal c to 
kagome planes. 6 g =36° denotes the tilt of CuCU pla- 
quettes with respect to the kagome planes [6], which 
locally sets the principal axes of the g-tensor. Typi- 
cal g-factor values for Cu 2+ ions in a uniaxial crystal 
field g\\ = 2.2 - 2.4 and g± = 2.05 - 2.15 should yield 
two distinctive narrow features in a single ESR deriva- 
tive spectrum (Fig.[TJ. In contrast, the observed spectra 
appear as a single broad derivative shape, which is a di- 
rect evidence of a sizable magnetic anisotropy. In the 
case of the DM interaction, which we argue below to 
be dominant, the symmetry of the lattice imposes the 
line-width to depend only on 0. To fit the spectra, we 
therefore convoluted a powder field distribution and a 
Lorentzian, the latter having the semi-phenomenological 
line-width AH = (AH Z cos 2 8 + AH p sin 2 6*) 1 / 2 , where 
AH Z and AH p correspond to the direction parallel and 
perpendicular to c. The field distribution combines a 
uniform angular distribution and the <?-factor anisotropy 
g{9). Standard g-tensor values, g\\ — 2.25 and g± = 2.14, 
were found. We stress that the RT ESR spectra are not 
noticeably affected by the Cu 2+ spins residing on the Zn 
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FIG. 2: (color online), (a) Scaling of ESR intensity Xesr and 
bulk susceptibility Xb in ZnzCu^z (OH)6Ch {x = 1). (b) The 
out-of-plane AH Z and the in-plane AH P ESR line-width. 



sites. They can be fitted with the same set of parameters 
for x = 0.5 — 1 (Fig. Q]). We therefore use them later to 
determine the magnetic anisotropy in the kagome planes. 

Several distinct regimes are observed in the In- 
dependence of the ESR line- width [Fig. EJJb)]. It 
is constant above ~200 K, with very similar values 
in ZnCu 3 (OH) 6 Cl 2 and Zn . 5 Cu 3 .5(OH) 6 Cl 2 . This In- 
dependence is characteristic of the high-T exchange 
narrowing regime. Below ^150 K a pronounced increase 
is observed in AH Z , while AH p starts increasing notice- 



ably below ^100 K. On general grounds [2l| . we recog- 
nize the increase of the line-width as the evidence of a 
building up of short-range spin correlations at T < J. 
The difference between AH Z and AH p indicates that the 
correlations are evolving in an anisotropic manner. We 
propose that this results from magnetic anisotropy af- 
fecting the evolution of the correlations. With decreasing 
temperature the ESR spectra should broaden monoton- 
ically as the spin correlations develop, e.g., similarly as 
found in another two dimensional (2D) frustrated system, 
SrCu2(BOs)2 [13|- At variance, a maximum is observed 
in AH Z . Its position shifts significantly toward higher 
temperatures with reducing the Zn content [Fig. [2(b)]. 
We therefore relate this maximum to the Zn 2+ /Cu 2+ 
inter-site disorder. The effect of different Cu 2+ spins on 
the ESR signal is proportional to their susceptibility [ijj. 
Therefore, below ~20 K where Xb overshadows Xk, the 
Cu 2+ spins on the Zn sites dominate the ESR response, 
which explains the crossover between 20 K and 150 K. 
Weakly coupled Cu 2+ spins on the Zn sites, with a differ- 
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ent environment, hence a different magnetic anisotropy, 
would provide a natural explanation of the low-T behav- 
ior. Substantial theoretical modeling is needed though 
to account quantitatively for the observed line-width In- 
dependence, which is beyond the scope of this paper. 

We address the origin and magnitude of the magnetic 
anisotropy in ZnCu3(OH)eCl2 from the detailed analysis 
of its RT spectrum, in the framework of the Hamiltonian 
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where, the first sum runs oven the nn spin pairs and rep- 
resents the exchange interaction H e , the second one gives 
Zeeman coupling Hz, and H' is magnetic anisotropy. For 
strong exchange with respect to Hz and H' the Kubo- 
Tomita (KT) formalism [23| yields a Lorentzian ESR 
spectrum, in agreement with our experiment. Regard- 
less of the specific form of the anisotropy, its line-width 
can be estimated from the magnitude A of the anisotropy 
as AH ~ A 2 1 J [3], yielding A ~ 16 K. 

We stress again that the spectra are extremely broad. 
Both, the dipolar interaction between Cu 2+ spins and 
their hyperfine coupling to nuclear spins are too small 
to account for the observed line-width. They would 
lead to AH < 1 mT due to strong exchange narrow- 
ing. Therefore, exchange anisotropy must be at work 
in ZnCu3(OH)eCl2. It originates from a sizable spin- 
orbit coupling A, which mixes orbital excited states of 
Cu 2+ with their Kramer's ground doublet, and is re- 
flected in measured g-factor shifts of 0.14-0.25 from the 
free electron value 2.0023. The second-order perturba- 
tion calculation yields two exchange anisotropy terms 
- antisymmetric DM interaction D oc (Ag/g)J (linear 
in A) and symmetric anisotropic exchange (AE) T oc 
{Ag/g) 2 J (quadratic in A) [IB]. The DM interaction 
H' = ^ij ■ Si x Sj is thus usually dominant in Cu- 

based antiferromagnets, if allowed by symmetry [25]. 

The out-of-plane DM component D z is generally 
present in the kagome lattice, while the in-plane com- 
ponent Dp is allowed in ZnCu3(OH) 6 Cl2 because the su- 
perexchange mediating 2 ~ ions break the mirror sym- 
metry of the kagome planes. We therefore base our ESR 
analysis on the dominance of the DM anisotropy. The 
Dy pattern (Fig. [I]) is obtained by choosing the direction 
of one of the vectors and applying symmetry operators of 
the lattice space group. We use the convention of spins 
being counted clockwise in all triangles. 

We can calculate the ESR line- width 1 26] 
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FIG. 3: (color online). The reduced x 2 , reflecting the quality 
of the RT fit in ZnCu3(OH)6Cb. The dashed rectangle gives 
the optimal parameters \D Z \ — 15(1) K and \D P \ = 2(5) K. 



determined by the DM pattern (Fig. [J), in the infinite-T 
limit. The high-T exchange narrowing regime is reached 
at RT since the ESR spectra do not change noticeably 
above ~200 K. This is consistent with a T-independent 
NMR spin-lattice relaxation above 150 K We fit- 

ted the RT ESR spectrum with the powder-averaged 
Lorentzian having the above line- width AH (9). In Fig. [3] 
we plot on the D z -D p map the reduced \ 2 obtained from 
the fit. Optimal parameters are \D Z \ — 15(1) K and 
| .Dp | = 2(5) K. According to Eq. <[2j) their sign cannot 
be determined by ESR. We find the magnitude of the 
extracted DM interaction (D/J w 0.08) in the range set 
by several other Cu-based 2D frustrated systems. For 
instance in the orthogonal-dimer system SrCu2_(B03)2 it 

Il,[27j|, while 



amounts to 4% of the isotropic exchange 
the ratio of 16% was reported for the triangular com- 
pound Cs 2 CuCl 4 [H]. 

At this point it is important to address the issue of the 
possible hidden symmetry of the DM interaction [29| . 
For linear spin chains with staggered DM vectors one 
can effectively transform it to the AE term of the mag- 
nitude D 2 1 J by applying a non-uniform spin rotation 
[3oT |. This makes it of the same order as the initial AE 
and discards the direct applicability of the KT formalism 
[3(1 [3l| . However, for other lattices one should distin- 
guish between reducible DM components, transforming 
to A 2 order, and irreducible components, remaining lin- 
ear in A [27] . The components that can be eliminated 
in the first order are those which sum up to zero within 
any closed loop on the lattice ji^. For the kagome lat- 
tice, the in-plane D p is reducible while the out-of-plane 
D z is irreducible. Using the KT formalism for the latter 
is therefore well justified, while the former might not be 
accurate. However, the small value of D v does not signif- 
icantly influence the ESR spectra (see Fig. [3]). We note 
that in SrCu2 (603)2 the values of the DM interaction 
obtained from the KT analysis [22], from inelastic neu- 
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tron scattering 27] and from NMR measurements [32| 
were found to agree within 20%. 

In Ref. [12] somewhat larger values, \D Z \ « 0.1J and 
\D p \ w 0.2 — 0.3 J, than extracted in this study were 
suggested to explain the bulk susceptibility upturn in 
ZnCu3(OH)eCl2, by considering the DM interaction as 
the only additional term to the Heisenberg Hamiltonian. 
These values would lead to high-T ESR line-widths of 
6—7 T, which is inconsistent with our experimental obser- 
vations. Therefore, DM interaction alone cannot account 
for the macroscopic susceptibility in ZnCu3(OH)gCl2. 

Sizable DM interaction (D/J ss 0.08) raises the fun- 
damental question to what extent DM anisotropy affects 
the low-T magnetism of ZnCu3(OH)gCl2. Spin-wave ex- 
citations in the spin-5/2 Fe-jarosite were explained by 
the presence of the DM interaction of a similar magni- 
tude (D/J « 0.09) [l3|. In this kagome compound, Neel 
order below 65 K was reproduced by classical Monte- 
Carlo calculations, predicting uniform (q = 0) ordered 
spin structures induced by the DM interaction [3]. In the 
5 = 1/2 case, quantum correction should suppress LRO 
for small D/J ||. For D/J > 0.05 recent numerical 
calculations seem to indicate a broken-symmetry state 
at T = 0, which however has no on-site magnetic mo- 
ment [33]. Further, the rather large Cu 2+ /Zn 2+ inter-site 
disorder may destabilize LRO. Indeed, vacancies on the 
spin-1/2 Heisenberg kagome lattice induce unexpectedly 
extended dimer-dimer correlations, despite short-ranged 
spin-spin correlations 0]. Although exact dimer patterns 
are not known for additional terms in the Hamiltonian 
that break the SU(2) symmetry, it seems plausible that 
they would contradict the uniform LRO state preferred 
by the DM interaction. This might explain the absence 
of LRO in ZnCu 3 (OH) 6 Cl 2 down to L»/300 = 50 mK. 

Last, quite importantly, the DM interaction mixes 
magnetic excitations into the presumably singlet ground 
state of the Heisenberg kagome lattice. This could 
explain the experimentally observed gapless excitation 
spectrum in ZnCu3(OH)eCl2 3, IB] and the finite sus- 
ceptibility of kagome planes at low T [13], as confirmed 
by numerical calculations on finite spin clusters [33 ]. 
These experimental results could likewise be described 
as a disorder-generated effect, due to a vacancy-induced 
spatially distributed density of states jjj^]. Most likely, 
both effects play an important role. 

In conclusion, we have shown that the DM interac- 
tion in ZnCu 3 (OH) 6 Cl 2 is sizable (D/J w 0.08) and 
could critically affect the low-T properties. Knowing 
its magnitude, future theoretical studies should investi- 
gate its impact on the ground state. They should also 
clarify whether the DM magnetic anisotropy and the 
Cu 2+ /Zn 2+ inter-site disorder are primarily responsible 
for the absence of the spin gap. 
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